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Abstract. We briefly review selected mathematical models that describe the dynamics of pattern 
formation phenomena in dip-coating and Langmuir-Blodgett transfer experiments, where solutions 
or suspensions are transferred onto a substrate producing patterned deposit layers with structure 
length from hundreds of nanometres to tens of micrometres. The models are presented with a fo¬ 
cus on their gradient dynamics formulations that clearly shows how the dynamics is governed by 
particular free energy functionals and facilitates the comparison of the models. In particular, we 
include a discussion of models based on long-wave hydrodynamics as well as of more phenomeno¬ 
logical models that focus on the pattern formation processes in such systems. The models and their 
relations are elucidated and examples of resulting patterns are discussed before we conclude with 
a discussion of implications of the gradient dynamics formulation and of some related open issues. 
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1. Introduction 

The deposition of complex fluids on solid substrates is a common part of various coating techniques 
employed for different purposes. Examples of such techniques are spin-coating, dip-coating, slot 
die coating and the doctor blade technique ll20l . Although the experimental processes are often 
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easily performed and also well controllable and reproducible, a detailed theoretical description is 
often still lacking, in particular, if complex fluids are deposited, where by complex we denote any 
fluid that has some ki nd of inner structure, e.g., a varying density of solute molecules. As a result 
pattern formation phenomena that commonly occur in such experiments Il29ll43l 68] are only rarely 
understood. This is due to the fact that among other issues a proper description of the deposition 
of a complex fluid on a substrate necessarily involves the modelling of the motion of a three phase 
contact line, i.e., of a region where substrate, complex fluid and vapour meet. For simple fluids, 
many different approaches and models describing such a contact line exist, however, it is still a 
field of active discussions and development [T6l f5f611. Accordingly, for a complex fluid a moving 
contact line poses additional problems and is therefore far from well understood. 

The importance of the contact line for pattern formation processes stems from the fact that 
structured depositions of a complex fluid or of one of its constituents normally occur in the prox¬ 
imity of a receding contact line which moves with respect to the substrate. Perhaps the simplest 
experimental setup to study this phenomenon is a drop of a complex fluid on a substrate, where the 
contact line moves due to forces resulting from capillarity and wettability and as well by evapora¬ 
tion/condensation 02],[29, 43[ 68]• This is, e.g., also the case for the well-known coffee stain effect 
[mmm 4S1 SOI • Although quite regular and also complex patterns can be achieved in such a set¬ 
up mmm, such a simple approach is, however, not very suitable for quantitative experiments, 
as important control parameters like the velocity of the contact line cannot be controlled directly 
but only indirectly (e. g., through a temperature-dependent evaporation rate). 

A simple way to externally influence (but not control) the contact line velocity is to employ a 
confinement that increases the (local) vapour pressure in the gas phase. An example is the so-called 
sphere on flat geometry, where the liquid is confined between a flat surface and a sphere on top of 
it m. In this geometry the liquid forms a capillary bridge between the substrate and the sphere. 
While the direction of the contact line movement is controlled in this set-up, the contact line speed 
is still only indirectly controlled via the influence on the evaporation process. 

A technique with immense practical impact is spin coating, where a drop of the fluid is placed 
at the centre of a substrate which is then rotated at high speeds. The acting centrifugal forces 
spread the fluid into a thin homogeneous film. The resulting film thickness can be controlled in a 
wide range via the rotation speed, however, the subsequent evaporation or dewetting process and 
therefore the contact line motion itself cannot be directly controlled in this set-up. 

Following the classification introduced in 11581 . all the aforementioned set-ups could be called 
passive , as there is no direct active control of the contact line motion. In contrast, e.g., the doctor 
blade technique introduces as an additional control parameter the velocity of the blade that removes 
excess fluid from a completely covered substrate at an adjustable thickness. For sufficiently small 
thicknesses of the resulting fluid layer, the contact line directly follows the velocity of the blade. 
A similar technique is proposed in II8T1 . where a polymer solution is deposited at the edge of a 
sliding glass plate above a resting plate, where both the velocity of the top plate and the distance 
between the two plates can be adjusted. This is also the case for slot die coating, where the liquid is 
deposited onto a substrate through a slot at an adjustable height above the substrate. Here the flux 
and the velocity with respect to the substrate are controlled. Note that most active set-ups can be 
used to produce either liquid layers with a small thickness that evaporate fast effectively coupling 


2 






M. Wilczek et al. 


Modelling Pattern Formation in Dip-Coating Experiments 


the imposed velocity and the contact line velocity or liquid layers with a large thickness. The 
latter evaporate slowly only after the actual coating process, making the contact line movement 
independent of the previous transfer velocity, effectively making it a passive set-up. 

In this work we focus on another active technique, the so-called dip-coating. As the name 
suggests, a substrate is coated by dipping it into a trough which is filled with the complex fluid 
that has to be applied to the substrate. The velocity at which the substrate is withdrawn from 
the trough has to be typically well regulated and represents an important control parameter of the 
experiment. The free surface of the fluid usually forms a meniscus where the substrate and the bath 
meet, which for small transfer velocities determines the properties of the contact line. Similarly to 
the doctor-blade technique or the slot die coating, dip-coating therefore provides a direct method 
to control the velocity of the receding contact line. Depending on the experimental conditions, 
various striped and branched patterns of different heights can be formed by this technique Il44l 45 ]. 
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Figure 1: AFM images of a DPPC monolayer on a mica substrate. The monolayer was transferred 
onto the substrate by Langmuir-Blodgett transfer with different transfer velocities v (direction from 
bottom to top) iflTTl . The different colours indicate different phases of the monolayer with different 
area densities and heights. 

A special case of dip-coating is the Langmuir-Blodgett (LB) transfer @|. For LB transfer, a 
trough filled with water is used, on which a floating surfactant layer is prepared. In this way a 
monolayer, i.e., a single layer of surfactant molecules, can be transferred onto the substrate. Due 
to the molecular interaction between the substrate and the surfactant film the surfactant layer can 
in the course of the transfer undergo a phase transition from a low density liquid expanded (LE) 
phase to a denser liquid condensed (LC) phase (64l[58l- This effect is called substrate-mediated 
condensation (SMC). Detailed experimental studies show that this transition can occur in a spa¬ 
tially inhomogeneous way, leading to various patterns formed by domains of the two different 
phases llb(T. ~28. 46. 111. The range of resulting patterns includes stripes of different orientations 
and lattice-like patterns. These patterns result from self-organisation processes as their type as well 
as their quantitative properties (e.g., the wavelength of stripe patterns) can be controlled through 
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experimental parameters related to out-of-equilibrium aspects as the transfer velocity of the surfac¬ 
tant layer or its density that is prepared at the water surface of the through. Note that, despite their 
striking regularity, these processes do not involve templates of any kind. Fig.[j]shows AFM images 
of a Dipalmitoylphosphatidylcholine (DPPC) monolayer on a mica substrate after LB transfer with 
different transfer velocities at the same surface pressure. 

With both, dip-coating and LB transfer, regular micrometre-scale patterns can be produced 
across macroscopic substrate regions. As this is a desired feature in a number of production pro¬ 
cesses, a thorough theoretical understanding of the occurring processes is needed. 

The structure formation in dip-coating experiments mainly occurs at the meniscus formed at the 
plate when it is withdrawn from the bath. However, to describe the physical processes in the region 
of the meniscus, one needs to formulate theories for the dynamics of thin films of different complex 
fluids on a solid substrate. In a computational fluid dynamics approach one would describe the 
system with a full-scale macroscale deterministic continuum model, namely, employing Navier- 
Stokes equations together with transport equations for the solute(s) and accounting for solvent 
evaporation and processes in the gas phase. This requests a large computational effort and suffers 
from uncertainties in the description of a number of processes. This approach is not reviewed here. 

However, for the considered geometries one is often able to develop asymptotic models - so- 
called thin film models. Here, the attribute thin does not necessarily refer to an absolute measure, 
but rather to a small ratio of typical length scales in directions perpendicular and parallel to the 
substrate Ii54l fT5l . Often the terms long-wave approximation and small gradient approximation 
are also employed in this context. For ultrathin films one may also apply alternative approaches 
based on microscale considerations as, e.g., ki netic Monte Carlo models or dynamical density 
functional theory lf72l [59 1. Here, however, we focus on the hydrodynamic thin film approach to 
the description of films of suspensions, mixtures and solutions. In particular, our brief review shall 
show that the various different systems of equations can be written in a unified way, namely, in 
a gradient dynamics form. This facilitates the comparison of the individual models and allows 
to better highlight their similarities and differences. It also allows to identify proposed models 
that do not comply with the gradient dynamics form and might therefore be thermodynamically 
inconsistent. 

Thin films of simple liquids (one phase, one component) are theoretically extensively studied 
and various long-wave models are well established [54] 15 , 66]. Several recent studies with such 
models consider the dip-coating geometry and conclude that even simple fluids show a very rich 
behaviour of the occurring meniscus structures |]62ll63lfT9, 82,,9l l26l 75]. However, the formulation 
of (long-wave) models for free surface films of complex fluids on solid substrates is less advanced 
and a present subject of intense research and indeed debate (see, e.g., [jT2] [13] 5], ESI, 25] 147 , 
m m no m n m mi). For a brief review of long-wave models see li68ll . Here we discuss 
the application of the underlying concepts to coating processes in active geometries J39] [38 , [40, 
46, 4F, 78J and briefly set them into the context of related approaches. The particular focus is 
the formulation of long-wave models in the unified gradient dynamics form that allows (i) for 
a straightforward check of thermodynamic consistency and (ii) for a systematic incorporation of 
additional physical effects. 
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First, we discuss in section [2~i~| the general form of gradient dynamics models for the evolution 
of a single field and the evolution of two coupled fields. Then we present two particular models 
for Langmuir-Blodgett transfer (section [2727]) and dip-coating (section [231] ) that can (i) be derived 
from the basic hydrodynamic transport equations and (ii) be brought into the general gradient dy¬ 


namics form for two coupled fields presented before. In addition, section 2.4. discusses a reduced 
model for the Langmuir-Blodgett transfer as well as possible extensions and generalisations. The 
work concludes in section [X] with a discussion of the relation of the presented models to related 
approaches. 


2. Mathematical Modelling 

2.1. General Form 

We begin with a short discussion of the general mathematical form of the conserved parts of the 
evolution equations, i.e., of the dynamics of a quantity that obeys a conservation law (or continuity 
equation). A well-known equation for a real-valued conserved order parameter u(x) is the so- 
called Cahn-Hilliard equation, that describes the dynamics of phase separation 0 driven by the 
particular underlying free energy functional T[u\ discussed in [[8|. In its general form, it reads 


d t u = -V • J = V • 


AT 

M {u V— 
ou 


( 2 . 1 ) 


with 




a -\Wu\ 2 + f(u) 


dx, 


( 2 . 2 ) 


where from now on we assume a two-dimensional domain flcK 2 and therefore x = (x. y) T e 
flei 2 and V = ( d x , d y ) v . although it can also be formulated for a domain of higher dimension. 
The free energy functional consists of a local free energy f(u) and an interface energy contribution 
11 Vu| 2 . As it is suggested by our notation, the Cahn-Hilliard equation has the form of a continuity 
equation, where the flux J = —M (u)'Vj^ is related to a thermodynamic force, represented by the 
gradient of a chemical potential — V|^, through the mobility function M(u). This ansatz for the 
flux is based on the assumption that the system is not too far from equilibrium. It can be easily 
seen (by calculation of t)] = f (SF/dufa t))d t u(x., t) dx), that the dynamics given by 

equation ( |2.1[ ) monotonically decreases the free energy functional T with time. 

Originally, the Cahn-Hilliard equation was derived as a model of nonequilibrium thermody¬ 
namics for the evolution of binary mixtures. However, equations of such gradient dynamics form 
appear in many places, e.g., they turn out to be of particular interest as hydrodynamic theories of 
liquid films on substrates in the thin film limit Il50l . There they can be derived from the Navier- 
Stokes equations with appropriate boundary conditions by employing the long-wave approxima¬ 
tion ll54l . The time evolution of the height profile //(x. t) of a thin liquid film on a substrate is 
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governed by the Cahn-Hilliard-type equation 


dth = V • 


3?/ Sh 


a , 


m= i — iv/if+ f(h) 


dx. 


(2.3) 


The cubic mobility function results when using no-slip boundary conditions at the substrate and 
stress-free conditions at the free surface and denoting the dynamic viscosity with q. In the free 
energy functional, the squared-gradient term represents the surface energy of the free surface of 
the liquid, where a is the liquid-gas interface tension. The local free energy f{h) is identified with 
a potential describing the interaction energy between the free surface of the film and the solid- 
liquid interface, which in the literature is referred to as disjoining potential, binding potential, 
adhesion potential or wetting potential, while the entire functional F[h] is sometimes referred to 
as effective interface Hamiltonian |[5]|. As the thin-film models in question do not contain terms 
due to inertia, i.e., they represent the overdamped limit, it is rather intuitive that the dynamics 
minimises a free energy functional through a Cahn-Hilliard type dynamics. Nevertheless, one has 
to bear in mind that although the free energy functional can be derived based on thermodynamic 
considerations, the mobility function can only be obtained from the long-wave expansion of the 
basic hydrodynamic equations. 

It is now only natural to formulate a similar model for thin films of complex fluids, which are 
described by n conserved order parameters u = {u\. u ->-.... «„) ' due to their inner structure. To 
describe the evolution of several order parameters governed by a single free energy functional, 


the Cahn-Hilliard equation ( |2.1[ ) can be generalised by replacing M (u) with an x n dimensional 
positive definite and symmetric mobility matrix Q(u) that encodes a thermodynamically linear 
coupling between all the thermodynamic forces — Vj^ and the fluxes of all order parameters. 
Introducing as well a symmetric nxn dimensional surface energy matrix £ in the free energy, the 
equations read 


d t u 


V • 


Q(U)V ^ 



1 

2 


(Vu) t -£(Vu) + /(u) 


dx. 


(2.4) 

(2.5) 


While u and J22,..., J22J are now //-dimensional vectors, one has to note that V and 

the dot product are still defined on the two spatial dimensions of position spacc0 The constraint 
2 The notation becomes clearer using an index notation with Einstein summation convention: 


d t u a = V 4 

^ [ul = / 
n 


Ca/3(u) V i 


51 " 
fiup _ 


2 (^/t/'ct) ^ct/3 


(ViU*) + /(u) 


dx, 


with * = 1,2 and a, /? = 1,.., n. 


( 2 . 6 ) 

( 2 . 7 ) 
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that the coupling of the thermodynamical forces — V|^ to all the fluxes has to occur through a 
positive definite and symmetric mobility matrix Q corresponds to the condition of positive entropy 
production and to the well-known Onsager reciprocity relations , respectively [52, 27]. 

As the systems that we discuss next are described by two order parameters, we explicitly give 


the expressions ( |2.4| ) and ( |2.5[ ) in the case n — 2: 
d t u = V • 

•Fful = 


Q(mi,m 2 )V — 

ou 


y |VUi| 2 + y |VU 2 | 2 + <7l 2 Vtli • V U 2 + f{u 1 ,U 2 ) 


dx. 


( 2 . 8 ) 

(2.9) 


An important mathematical consequence of the positive definiteness of Q is that the dynamics de¬ 


scribed by ( |2.8[ ) still monotonically decreases the free energy J~ over time. For the two-component 
case one has 


—J : [ui 1 U2] = 


57 dm 5T du 2 
5u\ dt 5m dt 


_ 6F\ _/< 5 T 


dx = [ ~^d t u dx 

J 5 u 

n 

dx < 0. 


( 2 . 10 ) 

( 2 . 11 ) 


n 


In the following sections, we present two models for structure formation during dip-coating which 
represent models for complex liquids of the mathematical form ( |2.8[ )-( |2l9| ). The presentation of 
the models in this unified gradient dynamics form allows for a direct comparison of the individual 
models and highlights their similarities and differences. Furthermore, one may argue that for over¬ 
damped (i.e., creeping flow) interface-dominated systems it should always be possible to provide 
such a gradient dynamics form that is automatically thermodynamically consistent. On this basis, 
on can then incorporate additional effects into the models in a thermodynamically consistent way 
by augmenting the free energy functional with additional contributions. 


Here, we treat two applications: section 2.2. considers the case of a liquid layer covered by an 


insoluble surfactant (see Fig. [2] (a) for a sketch and a definition of the variables), while section 2.3. 
treats films of a suspension/solution with a non-surface active solute (see Fig. [2] (b)). 


2.2. Langmuir-Blodgett Transfer 

As a first example for a two order parameter gradient-dynamics theory, we discuss a model for 
a thin film covered by an insoluble surfactant layer If69ll . This model can then be extended by 
advection and evaporation terms to model the transfer of a surfactant layer onto a substrate as it 
occurs in the Langmuir-Blodgett transfer described in the introduction |[39l . 

The natural order parameters to model a thin film covered by a surfactant are the density T(x) 
of the surfactant on the surface of the liquid film and the local height h(x) of the liquid film, see 
Fig. [2] (a). However, we need to identify order parameters which are conserved on the euclidean, 
non-curved surface A of the substrate. The physically intuitive F is a density with respect to a 
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(b) 



Figure 2: Sketch of the variables considered in the models for Langmuir-Blodgett transfer and 
dip-coating of a solution, (a): a pure liquid layer with the height /i(x) with a insoluble surfactant 
monolayer on top with a density T(x), as in the Langmuir-Blodgett transfer in Section 2.2. (b): a 
liquid solution layer of the height /i(x) with a solute density <^(x), as in the dip-coating system in 
Section 1231 


curved film surface element dS, i.e., Y and h cannot be varied independently. A conserved order 
parameter pendant to F is the density 


f = £r, with f = \J 1 + ( d x h ) 2 + ( d y h ) 2 , 


( 2 . 12 ) 


which is the projection of the density field on the curved film surface onto the surface of the 
substrate, i.e., dS = £dx. To model a system, where a surfactant phase transition may occur, 
we employ a local free energy that depends on Y and may, e.g., allow for two coexisting phases. 
From this local free energy follows the dependence of the surface tension of the film on surfactant 
density. A general free energy functional for such a system in terms of the independent conserved 
order parameters reads 


T[h,Y] 


f{h) + g 






(2.13) 


In this functional, f(h) is the wetting energy of the film as for the case of a simple liquid in 
Eq. ( |2.3[ ). The function g(Y, h) is the local free energy of the surfactant layer, which allows for a 
phase transition which can be triggered by a change in film height. We use 


g{ r, h) = Y, G k (r - r cr ) fc + Bf(h ) (r - r cr ), (2.14) 

k=0 

where the Gk and B are constants. The first term is a double well potential, i.e., a fourth order 
polynomial in Y centred about a critical density r cr - a standard model for a first order phase 
transition. Note that the zeroth order term Go corresponds to the surface energy of the liquid layer 
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with a homogeneous surfactant coverage of density r cr (denoted as a in Eq. ( ]2.3[ )). The final term 
in g(T, h) is linear in (T — r cr ) and models the influence of the substrate-mediated condensation 
(SMC) in the form of a film-height-dependent tilt of the double well potential, that energetically 
favours the condensed phase when the distance between the surfactant and the substrate, i.e., the 
height of the liquid film, becomes small. There exist different ways to implement the tilt, in 
particular, as no measurements of surfactant pressure isotherms in dependence of film heights are 
available. It is known, however, that the SMC lowers the free energy of the high-density LC phase 
for a vanishing distance between the monolayer and the substrate. Furthermore, this influence 
should decay fast with increasing film heights and disappear over a length scale comparable to 
the length scale over which the disjoining potential is acting. Therefore, here the latter’s height 
dependence is chosen to control the tilt ||39|| , but other choices are conceivable. 


The final term in Eq. ( |2.13| ) is an energy contribution for the surfactant that penalises spatial 
inhomogeneities of the surfactant density, i.e., it is the energy of the one-dimensional interfaces 
between areas on the film surface with different surfactant phases. The parameter k is a corre¬ 
sponding interface stiffness. Note that the various factors £ in the free energy functional T stem 
from the projection of surface densities and surface gradients onto the substrate plane. 

The full gradient dynamics is given by 


d t h = V • 
d t f = V • 


„ 

Qhh V — + Q hr V —~- 

Oh ST _ 

Qvhy tt“ + Qrr V —~g 

Oh 5T 


(2.15) 


where the mobility matrix is 



(2.16) 


with the molecular mobility D. Note, that although the mobility matrix has always to be obtained 
from hydrodynamics, one may generalise the set of equations by amending or extending the free 


energy functional. The equations one obtains when performing the variations of (2.13) are (i) in the 


low concentration limit (i.e., using Eq. (2.13) with entropic contributions to g(T) only, and with 
k = 0) f[69ll identical to equations one obtains via a long-wave approximation of an advection- 
diffusion equation for the surfactant and the Navier-Stokes equations with appropriate boundary 
conditions at the substrate and at the surfactant covered surface (where a linear Marangoni effect is 
incorporated) [35, J 5] and (ii) also recover the model for substrate-mediated condensation derived 
in 


There, the equations resulting from ( |2.15[ ) with the free energy functional (2.13) were 
directly derived based on a long-wave approximation and then supplemented with an evaporation 
term and an advection term to model the Langmuir-Blodgett transfer]^] In order to incorporate the 


3 The resulting equations contain a term d/,g(l'. h) that was missed in the derivation presented in 
leads to a minor rescaling of the wetting energy. See Ref. fl69l for details. 
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evaporation term in a fully consistent way with the free energy functional, the term 


5 


ev 


—E„ 



(2.17) 


may be added to the evolution equation for h. This corresponds to an extension of the conserved 
Cahn-Hilliard type equation by a non-conserved Allen-Cahn type contribution. The extension 
is based on the assumption that the evaporation is a close to equilibrium process driven by the 
difference of the chemical potential of the solvent layer given by and the chemical potential fi v 
of the vapour phase which is assumed to be constant and homogeneous. The parameter E v is a 
constant that defines the time scale of the evaporation compared to the time scale of the conserved 
dynamics. The full resulting equations then read (cf. 1139116911 ). 


d t h = V • 
d t T = V • 


[d h f(h) - V ■ ( 7 Vft)] - ^V 7 - h\ 


T ^ev) 


h 2 Y 

2r] 


(hr 


V [d h f{h) - V • (7V/1)] — ( — + D ) V7 — Tv 


(2.18) 

(2.19) 


where we have introduced a generalised ’surface tension ’7 = g - r<9 r 9 - f (vr ) 2 + kTV 2 t, 
as well as the transfer velocity v. Note, that in Il39l the approximation V ■ (7V/1) « 7V 2 /! 
was used. Direct numerical simulations of the model completed with suitable lateral boundary 
conditions (see [f39ll for details) show good qualitative agreement with the experimental results. 
The range of patterns that can be obtained is illustrated in Fig. [3] We emphasise that the general 
formulation reviewed here allows for a straightforward incorporation of additional physical effects 
in a thermodynamically consistent manner and to identify proposed models with ad-hoc additions, 
e.g., for surfactant-concentration dependent disjoining pressures, that lack this consistency. Several 
examples are discussed in [f69ll . A further example is an additional control of the pattern formation 
that can be achieved by the use of prestructured substrates. This can be modelled by a spatially 
modulated wetting energy f(h, x) and results in the occurrence of more complex patterns including 
locking effects between periodic patterns and the prestructure on the substrate [|38l . 
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Figure 3: Snapshots from direct numerical simulations of the model for Langmuir-Blodgett transfer 
( |2.15[ ) [39j [37’], where the height indicates the height of the liquid and the colour scheme encodes 
the density of the surfactant layer. Depending on the transfer velocity v, one observes the depo¬ 
sition of a homogeneous high-density layer, a stripe pattern orthogonal to the meniscus, a stripe 
pattern parallel to the meniscus, or a homogeneous low-density layer. For simulation details and 
parameters see ll37l . 

2.3. Dip-coating 

As second example for a two order parameter gradient-dynamics model, we next discuss the case 
of a film of a solution of molecules or of a suspension of (nano)particles (that are assumed to be not 
surface active) that is transferred onto a substrate via dip-coating. Similar to the case of LB transfer 
described in the previous section, a thin film model for dilute solutions may be derived starting from 
hydrodynamic equations (Navier-Stokes, diffusion-advection equation). At first sight, it might 
seem that the natural order parameters for such a model are the height /r(x, t) of the film and a 
height-averaged concentration of the solute, 


h(x,t) 

^( X,t ) = h(x t ) / c ( x ’^) d ^ ( 2 ‘ 2 °) 

o 

where x = ( x , y) are again Cartesian coordinates defined on the surface of the substrate and 
c(x, z, t ) is the standard bulk concentration. A thin film model for the time evolution of h and 
4> was derived in li76l and also used, e.g., in [72J, 23, [24], while a gradient dynamics form in 
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terms of variations in h and 0 is given in l[T3Tl . However, the fields h and 0 should not be varied 
independently as for a fixed local amount of solute a variation in height results as well in a change 
in 0. Ref. lH3ll accounts for this by employing constrained variations (for a brief discussions of 
problems we see with this approach see 1167116810 . 

A simpler approach is taken in [67, [Tl), where it is argued that the natural independent order 
parameter fields are the film height /i(x, t) and the local amount of solute T(x. t) = //(x. t)0(x, t), 
i.e., an effective local solute ’layer height’. Based on /i(x, t ) and T(x, t ) the dynamics may be for¬ 
mulated as a gradient dynamics based on a free energy functional that corresponds to an extended 
interface Hamiltonian. In terms of the conserved order parameter fields h and T. the equations for 
a thin film of solution/suspension read [67I17T| 



SJ 7 SJ 

d t h = V • 

«mV s +«„V s 


5J 7 „ AT 

•e 

<1 

Q * kV jh +QnV m 


with the mobility matrix 


Q = 


h? h 2 V 

31 ; 3?7 

h0. ^ + m 


and the free energy functional 




/ 


f(h) + ^\Vh\ 2 + hf B 



dx. 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


Here, the first two terms are the respective wetting and surface energy as for the simple liquid 
film. The third term corresponds to the bulk free energy of the solute integrated over film height, 
i.e., it contains the concentration dependence of the free energy. In the dilute limit (i.e., including 
only entropic contributions valid at low solute concentration, ftiio) ~ 01n(0)) one recovers the 
standard hydrodynamic thin film equations that may then be expanded systematically by incorpo¬ 
rating further entropic terms and interaction terms into the free energy functional, analogous to 
the classical approach of Cahn and Hilliard [|8l, or simply using the Flory free energy for polymer 
solutions If22l . One may also incorporate influences of the solute on wettability [67i 17Tl . A deriva¬ 
tion of a fully solute-solvent symmetric version on a route via Onsager’s variational principle ll53ll 
is discussed in [f80l . 

As an example we consider the incorporation of entropic effects over the entire concentration 
range for the solute and include a Flory interactions term. Then the bulk free energy reads 

fstt 0 = ^0 [01n(0) + (1 - 0)ln(l - <#>)] + cap (1 - f) + ||V<#>| 2 = m + - 2 \V0\\ (2.24) 

where a is a microscopic length scale related to the solute and a is the Flory interaction parameter. 
The equations resulting from such a free energy for a general g with the additional gradient term 
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Figure 4: Snapshots from direct numerical simulations of the model for dip-coating ( |2.21[ ), where 
the height indicates the height of the liquid and the colour scheme encodes the height-averaged 
concentration 0 of the solute. Depending on the initial concentration 0 O , the transfer of drops 
in a hexagonal pattern or stripes parallel to the meniscus is observed. For simulation details and 
parameters see 
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in 0(x, t) are discussed in IfTPl . The specific form in Eq. (2.24) is discussed in Il65ll . To model 
the transfer process in dip-coating, we augment the gradient dynamics equations for conserved 
fields (2.21) with an appropriate evaporation term and an advection term to model the transfer 
due to dip-coating as done before in section [2^2] In particular, evaporation results in an increased 


concentration within the front region of the meniscus, but should also saturate for concentrations 
approaching unity. Following the same arguments as in section [2T2T[ we choose for the evaporative 
contribution 


- l-ivj ■ (2.25) 

Evaluating the variation of J~ with respect to h, while keeping T fixed, gives g — og'. This is a 
contribution of osmotic pressure type to the chemical potential of the solvent. Note, that without 
the constraint of fixed solute amount only g' would result. Osmotic pressure type contributions to 
evaporation are not frequent in the literature and are often introduced in an ad-hoc manner lf68l . 
However, in the gradient dynamics formalism it emerges naturally. For instance, when the entropic 
contribution ~ (1 — 0)ln(l — </>) is incorporated into g, the osmotic pressure contribution ensures 
that evaporation saturates for </>—>■ 1. Including advection terms and appropriate lateral boundary 
conditions, one may study the pattern formation occurring during dip-coating experiments [[651 . 
Direct numerical simulations of the presented model show various self-organised pattern types, 
such as, stripes parallel to the meniscus and hexagonal patterns, whose details and transitions 
depend on the transfer velocity and the initial concentration of the solute. Figure [4] shows two 
snapshots from simulations with different initial concentrations, exhibiting the two mentioned pat¬ 
tern types. Similar patterns were also found in Iflffl . where a related active transfer geometry is 
experimentally investigated. 


S„, = -E„ ( ^ 


2.4. Reduced Models 


Under some conditions it is justified to simplify the models presented in the previous two sections. 
Besides eliminating contributions to the free energy that account for a physical effect or interaction 
that turns out to only have a minor impact in the experimental system, one may also try to reduce 
the number of order parameter fields describing the system. The latter approach is feasible in 
the case of the model for the Fangmuir-Blodgett transfer (|2.15[). Time simulations of the model 


indicate that the film height fi(x, t) relaxes on a rather short time scale onto an almost stationary 
state, whose exact form has no major influence on the surfactant dynamics. Therefore it seems 
reasonable to replace the dynamic equation for the film height with a static approximation that 
enters the equation for the surfactant density T only parametrically. In addition, one can show that 
a static height profile minimises the generalised pressure gradient, making it possible to neglect 
the corresponding term as well in the evolution equation for the surfactant density ll40]| . The 
resulting equation for the surfactant density is of simple Cahn-Hilliard type (see section [TF]) with 


an additional advective term. It can be written in the gradient dynamics form ( |2.1[ ), where we now 
use c(x, t ) as the order parameter for the surfactant density to avoid confusion with the full model 
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(2.15): 


d t c = V ■ 


M(c)V- 

oc 


v.Vc, F[c] = J 


2 lV<f+ /(c) 


dx, 


(2.26) 


with the mobility M (c) and the local free energy /(c) given by 

M{c) = 1 and /(c) = ^c 4 - ^c 2 + /<(x)c, 


(2.27) 


respectively. The local free energy /(c) is modelled as a double well potential with two local 
minima corresponding to the two phases of the surfactant layer, similar to ( |2. 14[ ). The linear term 
introduces a tilt of the potential favouring the high concentration phase that is proportional to a 
spatially inhomogeneous term ('(x), which parametrically mimics the shape of the static liquid 
meniscus occurring in the full model [[40]]. The resulting reduced model is considerably less com¬ 
plex, but still captures the essential physical aspects of the full model. This can easily be seen by 
direct numerical simulations of the model, see Fig. [5} The resulting patterns are intriguingly similar 
to the ones obtained with the full model, which in turn clarifies the origin of the pattern formation 
process to be the phase separation dynamics of the surfactant monolayer, while the hydrodynamics 
of the underlying liquid layer is secondary. 



Figure 5: Different types of patterns obtained by direct numerical simulation of the reduced Cahn- 
Hilliard model for Langmuir-Blodgett transfer ( |2.26[ ). The range of patterns comprises stripe pat¬ 
terns with different orientation and adjustable wavelengths Il40l ]78 |. For simulation details and 
parameters see 


In addition to direct numerical simulations, the simplicity of the reduced model allows for 
a more thorough analysis. A detailed bifurcation analysis of stationary states of the model was 
presented in ll40ll . also identifying the underlying local and global bifurcations that trigger the 
pattern formation process. The emergence of the entire bifurcation structure itself was furthermore 
investigated in [|4T1 . 

The reduced model lends itself to the introduction of additional effects, e.g., an interaction of 
the pattern deposition with a substrate prestructuring. As in the case of the two-component model 
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Figure 6: Snapshot of a direct numerical simulation of the reduced Cahn-Hilliard model for rotat¬ 
ing Langmuir-Blodgett transfer ( |2.26| ) modelled by using an advection velocity field v(x, y) that 
corresponds to a rotation of the substrate around a centre at ( x c , y c ) = (0, —500) with a rotation 
frequency of oj = 4.6 • 10 -5 , leading to different effective transfer velocities across the substrate. 
As a result, the pattern type varies across the substrate from (i) a homogeneous high density area to 
(ii) stripes perpendicular to the meniscus to (iii) stripes parallel to the meniscus to (iv) a homoge¬ 
neous low density area If77l . Due to the strong velocity gradients, the stripe patterns tend to break 
up into droplets aligned on a line. For further simulation details and parameters see [1771 . 


for LB transfer ||38l . such an interaction can lead to the production of more complex pattern types 
and to the occurrence of locking effects between the periodic prestructure and the resulting patterns 

m. 

Besides modifications of the free energy, also changes of the transfer geometry and process can 
easily be incorporated, e.g., one may consider the so-called rotating transfer, where the substrate 
is not drawn out of the trough in a straight way, but instead with a rotational movement [TTOl . With 
this technique, different parts of the substrate are withdrawn at different transfer velocities, leading 
to a non-constant v(x). As the pattern formation process is strongly influenced by the transfer 
velocity, this yields complex patterns such as the one shown in Fig. [6] which can be understood as 
a sequence of the basic pattern types presented in Fig. [5]that all exist on a single substrate [|77l . 


3. Conclusion 

We have briefly reviewed three models that describe pattern formation phenomena in experiments, 
where complex fluids are transferred onto a substrate producing patterned ultrathin deposit layers 
with structure length from hundreds of nanometres to tens of micrometres. All models have been 
presented within the framework of a gradient dynamics formulation that clarifies how the various 
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dynamics are governed by particular free energy functionals. These may be seen as extensions 
of the interface Hamiltonians for simple liquids. While the presented two-field models can be 
derived from basic hydrodynamic equations using a thin film approximation, we emphasise that 
the re-formulation as a gradient dynamics enables one on the one hand to easily compare and relate 
the different models to each other. On the other hand it allows for a direct and consistent way to 
extend the models by amending the individual contributions of the free energy and by including 
additional ones that account for further interactions and effects. It also offers a way to classify all 
such long-wave models derived from hydrodynamics based on whether they can be formulated as 
gradient dynamics or not. 

Finally, we would like to highlight some open issues and relate the presented models to se¬ 
lected literature results in the wider field. Inspecting Eqs. ( |2.15[ ) with ( |2. 1 6[ ) and Eqs. ( |2.2 1 [ ) with 
( |2.22[ ), one might even see these formulations as ’trivial’ because they have exactly the form that 
is expected for the evolution of two scalar conserved order parameter fields in the context of linear 
nonequilibrium thermodynamics Il52ll53ll . However, for two-field long-wave hydrodynamic equa¬ 
tions this had to our knowledge only been shown before for dewetting two-layer films [55] 156] 
but neither for surfactant-covered films nor for films of solutions or suspensions. Recently, the 
gradient dynamics form in the case of two-layer films has been extended to include the case with 
slip and has been employed in deeper mathematical analyses of such models [32j [34, 33J. We 
expect that the here reviewed gradient dynamics form will also facilitate similar analyses for films 
of suspensions, solutions and mixtures. 

Note that the request that similar models for films of complex fluids in relaxational situations 
(without additional external driving forces and without inertia effects) should have a similar gra¬ 
dient dynamics form provides one with a good test for model consistence that a number of ad-hoc 
amendments introduced in literature models do not pass. It also provides one with a guide how 
to proceed in the further development of multi-field models, e.g., it is an interesting question how 
long-wave models for surface-active solutes and soluble surfactants [[76l IT5 1 can be brought into 
gradient dynamics form even in the low concentration limit. The next step could then be to intro¬ 
duce into the relaxational gradient dynamics models additional nonequilibrium driving forces that 
break the gradient dynamics form in a controlled way - an example are liquid films covered by a 
carpet of self-propelled surfactant particles [fll [571. 

Another interesting aspect is that the presented gradient dynamics form brings the long-wave 
hydrodynamic models into the context of models known as dynamical density functional theories 
(DDFT) for the time evolution of the density field(s) of (colloidal) fluids. They represent (normally 
diffusive) gradient dynamics on equilibrium free energy functionals derived with the methods of 
statistical physics ll48l l2ll. i.e., their main underlying assumption is that the two-point correlation 
function out of equilibrium can be represented by its equilibrium variant ll48ll . This implies that 
the evolution equation for one and two fields is of the form of the present Eq. (2.1) and Eqs. ( |2.8 ), 
respectively. However, due to the character of the energy functionals, the evolution equations are 
often integro-differential equations instead of the partial differential equations encountered here. 
Note that in the two-field DDFT introduced in lf72l [591. the diffusive mobility matrix is normally 
diagonal. How this is related to the choice of frame of reference for diffusion is discussed in the 
final part of lf80l . It will be interesting to further pursue consequences of the formal similarity of 
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DDFT and thin-film hydrodynamics, e.g., to develop models that combine the various pathways of 
transport. 
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